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Abstract 

In this work we have analysed some cosmological bounds concerning an open FLRW solution of 
massive gravity. The constraints with recent observational H (2) data were found and the best fit 
values for the cosmological parameters are in agreement with the ACDM model, and also point 
to a nearly open spatial curvature, as expected from the model. The graviton mass dependence 
with the constant parameters 03 and <24, related to the additional lagrangians terms of the model, 
are also analysed, and we have obtained a strong dependence with such parameters, although the 
condition rn g ~ HJ 1 seems dominant for a long range of the parameters 0:3 and 0:4. 
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I. INTRODUCTION 


Current observations of Supernovae type la (SNIa) [11,12], Cosmic Microwave Background 
(CMB) radiation @, and Hubble parameter data pfl indicate an accelerated expansion 
of the universe, being the ACDM model the best model to fit the observational data. The 
A term corresponds to a cosmological constant energy density which is plagued with several 
fundamental issues which has motivated the search for alternatives models of gravity 
that could explain the observations. 


Massive gravity theories [s 


20 ] are old candidates to explain the accelerated expansion 


of the universe, since that the graviton mass could perfectly induce and mimic a cosmo¬ 
logical constant term. However, such kinds of theories were considered for long time as 
being unsuitable due to the appearing of Boulware-Deser (BD) ghosts js]. Recently it was 
discovered a nonlinear massive gravity theory that was first shown to be BD ghost free by 
Hassan et al. |10,|ll|, and also in the Stuckelberg formulation in [12]. Then such theory was 
also developed by de Rham et al. [13], sometimes called dRGT (de Rham-Gabadadze-Tolley) 
model (see 14] for a review), bringing back the cosmological interest for such theories 15]. 
Self-accelerating cosmologies with ghost-free massive gravitons have been studied thereafter 
US i- Nowadays there is a general agreement that in dRGT models which are defined 


with a flat reference metric, isotropic flat and closed FLRW cosmologies do not exist, even 
at the background level. Nevertheless isotropic open cosmologies exist as classical solutions 
but have unstable perturbations. In the case of a non-flat reference metric a ghost free 
theory exist for all types of background cosmologies ll|, though the perturbations are still 
unstable in isotropic cases [19]. 

Although the theoretical aspects of massive gravity has been severely studied in the 


last years, the cosmological constraints with observational data does not. In 24] has been 


investigated the cosmological behavior in the quasi-Dilaton nonlinear massive gravity, and 
the parameters of the theory has been constrained with observational data from SNIa, 
Baryon Acoustic Oscillations (BAO) and CMB. 

In order to study some cosmological bounds concerning the free parameters of the theory, 
we must have a massive gravity theory that has a FLRW limit well established for the 
reference metric. As showed in 19] the perturbations in unisotropic FLRW are stable, thus 
if we suppose that the unisotropies that cure the instability of the model do not greatly 
affect the cosmic evolution, we can use the results involving isotropic solutions as a first 
approximation to the stable unisotropic case. For this case the results of pj| are a good 
starting point to study some bounds in the massive gravity theory compared to recent 
observational data. In Section II we present the general massive gravity theory and the 
cosmological equations in Section III. In Section IV we present the constraints from H(z) 
data and in Section V some bounds on the graviton mass are presented. We conclude in 
Section VI. 
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II. MASSIVE GRAVITY THEORY 


Our starting point for the cosmological analysis is the massive theory for gravity proposed 
in 13]. The nonlinear action besides a functional of the physical metric g, iv {x) include 
four spurious scalar fields (p a (x ) with a = 0,1, 2, 3 called the Stiickelberg fields. They are 
introduced in order to make the action manifestly invariant under diffeomorphism, see for 
example 2l|. Let us start by observing that these scalar fields are related to the physical 
metric and enter into the action as follows: 


9i±v — ffiu T H^ u (1) 

where it is defined the fiducial metric f, w which is written in terms of the Stiickelberg fields: 

Uu = fa b md,rd^ b . ( 2 ) 

Usually f a b is called the reference metric and for the purpose of this paper, one can use 
fab = Vab = (—, +5 +5 +) once the scheme proposed by dRGT respects Poincare symmetry. 
Hence the fiducial metric in ([I]) is nothing but the Minkowski metric in the coordinate 
system defined by the Stiickelberg fields. So we have automatically defined the covariant 
tensor H which propagates on Minkowski space and the action is then a functional of the 
fiducial metric and the physical metric g^ v . 

The covariant action for massive general relativity that we are going to work with, can 
be written as: 


5 = M 2 Pl 


J d 4 x \f^g 


^ + 7n 2 U(g, H) 


( 3 ) 


where U is a potential without derivatives in the interaction terms between H^ v and g, liy 
that gives mass to the spin -2 mode described by the Einstein-Hilbert term. As observed 


m 


13] a necessary condition for the theory (J3]) to be free of the Bouware-Deser ghost in 


the decoupling limit is that \J—g U(g,H) be a total derivative. The most general covariant 
mass term which respects this condition is composed by 1 : 


J d A x yj^g U(g , H) = j d^xyf^g {a 2 C 2 + a 3 £ 3 + o^Cf) 
where a, are constants and the three lagrangians in (J4]) are written as: 


( 4 ) 


C 2 = -([/C] 2 -[/C 2 ]) 


( 5 ) 


20( a most general formulation is presented, where the parameters a, are assumed to be dependent on 


In 


the Stiickelberg fields. 
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( 6 ) 

( 7 ) 


C 3 = i([/C] 3 -3[/C][K 2 ] + 2[K s ]) 

= _L([/c] 4 -6(K] 2 [/c 2 ] + 3[x: 2 ] 2 + 8[x:][K 3 ]-6[K 4 ]), 


where one has defined the tensor /C(f = 8^ — \/d^(j) a d u (j) b pab- In general there are other 
polynomial terms in /C and the procedure to generate it can be found in fill. 131, however 
it has been shown that all terms after the quartic order vanishes (see also 22 


,2 


Notice 


also that we have maintained the constant term 02 in order to see its consequences on the 
evolution equations. The case 02 = 0 is pathological since in this case the linearised theory 
and the non-linear one have different number of propagating modes. 


III. COSMOLOGY OF MASSIVE GRAVITY 


Let us begin by considering an open (K < 0), homogeneous and isotropic FRW universe 
for the physical metric: 


g^dx^dx u 


— N(t) 2 dt 2 + a(t) 2 


dxidx 1 + 


K{xidx 1 ) 2 
1 — IiXiX i 


( 8 ) 


where, p, u = 0,1, 2, 3 and i, j = 1, 2, 3, with a; 0 = t, x 1 = x, x 2 = y , x 3 
same ansatz for the Stiickelberg fields used in 19], i.e: 


z. Adopting the 


'. 7 1 = f(t)V 1 - K ; O' = V-Kf(t)x\ 


( 9 ) 


after plugging back the metric and ([9]) in ([3]), one obtains the following Lagrangian 
a{t ) and f(t), where overdot will denote the time derivative: 


19] for 


C 


9 


1 

8nG 


3KN(t)a{t) 


3d (t) 2 a(t) 

W) 


+ 771^(0:2^2 + 03L3 + 04L4) , 


( 10 ) 


where 


L 2 = (3a(f) 2 - - f(t)a(t) - N(t)V=Kf(t )), 

L 3 = (a(t) - V^Kf(t)) 2 (W(t)a(t) - 3 f(t)a(t) - Ny/=Kf{t )), 

L 4 = ( a {t) - V—Kf (f)) 3 (V (t) - /(f)). (11) 

The matter content is assumed to be of the form Tjf = diag[— p m (t), p m (t), Pm(t), p m (f)]. 
Taking the Euler-Lagrange equation of C g with respect to / leads to 

(d(t) - y/^KN(t))[ o 2 (3 - 2 C) + o 3 (3 - 4C + C 2 ) + o 4 (l -2C + C 2 )} = 0 , (12) 
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where C = f(t)y/—K/a(t). The two interesting solutions are 2 : 


C± = 


o?2 + 20:3 + 0:4 A a/AJ 

(0:3 + 0:4) 


( 13 ) 


where A = 0 : 2(02 + 03 — 04) + a 2 . We take Euler-Lagrange equation of (fTOf) with respect to 
N and use (1T3|) to obtain the Friedmann equation: 

a(t) 2 K 87 tG A± 

N(t) 2 a(t ) 2 a(t) 2 3 Pm 3 ’ 1 ’ 

where A± = m 2 g l3± and 


At 


1 

(03 + 04) 2 


202+302(03—04±\/A)+3 o 2 (o 3 —o 4 )(o 3 ±\/A)—(o 3 ±\/A) 2 (— 2 o 3 ±v^A) 


(15) 

is a dimensionless parameter depending only on the constants 02, 03 and 04. In this equation 
one can recognize A± as being the energy density of massive gravity (p g ). Equally, it is 
expected to find the pressure term (p g ) in the second Friedmann equation. 

Such equation can be obtained by combining Eq. (TT4|) with the variation of Eq. flTOjl with 
respect to a{t). Namely, we subtract Eq. (TT4|) from the Euler-Lagrange equation of a{t) in 
order to get 


2 H{t) 2 K _ 87 tG ( 


(16) 


where H = • Notice that the r.h.s of the above equation contains only the contribu¬ 

tion from the matter part, which indicates that the graviton mass contribution satisfies an 
equation of state of the form p g = —p g , exactly as a vacuum behaviour. This same result 
was observed in (l9| . 

Another combination is possible in order to get a direct relation among a(t), pressure and 
energy density of matter and graviton. At this time we eliminate w of the expression 

from variation of Eq. (flOl) with respect to a(t). Thus, it is possible to get 


a(t ) 

N(t)a(t) 


A± 

3 


4v tG 

~1T 


(Pm + Ap m ). 


(17) 


With such equation it is much easier to analyse the universe acceleration. We conclude that 
an accelerated expansion occurs when A± > 47rG(p m + 3 p m ). 

It is also easy to see that A± acts exactly like an effective cosmological constant in Eq. 
(HU). In both Friedmann equations, (fill) , (flfil) and also (1171) . there should be set N = 1 in 
order to reproduce a cosmological scenario. In order to reproduce a positive cosmological 


2 The case a(t) = \/ — KN(t) just reproduces a constant scale factor when we take N —> 1 in order to recover 
the FLRW metric. 
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constant (which leads to an accelerating universe), we must have A± > 0, which implies 

£±>o. 

From now on we will assume a 2 = 1 according to the original dRGT theory 13]. The 
Friedmann equation (1141) can be rewritten in terms of the present critical energy density 
p c = 3 Hq/8itG, 

H {t) 2 = H o (—) + ™, 2 P± - \ , (18) 

V p c J y cr 

where H{t) = a/a is the Hubble parameter and H 0 ~ 70 km/s/Mpc is its present day value. 
Writing p m = p m o( a /a 0 )~ 3 , where p m0 is the present day value for the matter energy density, 
and introducing the density parameters 


fim — 


PmO 

Pc 


fig = At 


m 9 

H 2 ’ 


n K = 


K 


, 2/72 ’ 

i 0 -n 0 


the Friedmann equation (fT8]l can be expressed as 


Hitf = H { 


n r 


aoY 
a J 


+ fi K 


aoY 
a J 


+ fia 


or in terms of the redshift parameter, defined by 1 + z = a^/a, 

H{z) 2 = Hq fi m (l + z) 3 + (1 — fi m — fi 9 )(l + z) 2 + fi s 
where we have used the Friedmann constraint 


(19) 


( 20 ) 


( 21 ) 


1 — + fi/f + fi g , (22) 

that follows from (fT8jl and ffl9]l . Observational data can be used to constrain the values of 
such parameters and this will be done in the next section. 


IV. CONSTRAINTS FROM OBSERVATIONAL H(z) DATA 


Observational H{z) data provide one of the most straightforward and model independent 
tests of cosmological models, as H{z) data estimation relies on astrophysical rather than cos¬ 
mological assumptions. In this work, we use the data compilation of H(z) from Sharov and 
Vorontsova Q, which is, currently, the most complete compilation, with 34 measurements. 

From these data, we perform a y 2 -statistics, generating the Xh function of free parame¬ 
ters: 

34 

Xh = X] 

i =1 

where E{z) = and H{z) is obtained by Eq. (j2T]h 

As the function to be fitted, H(z) = H 0 E(z), is linear on the Hubble constant, H 0 , we 
may analytically project over H 0 , yielding x 2 H : 


H 0 E( Zi , fi mi fig) Hi 


v Hi 


(23) 


Xh = C~ 


B 2 

~A 


(24) 
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FIG. 1: Solid lines: Statistical confidence contours of massive gravity from H(z) data. The regions 
correspond to 68.3%, 95.4% and 99.7% c.l. Dashed line: Flatness limit, where fl m +fl g = 1. Points 
above this line are not considered on the statistical analysis. Star point: Best fit, corresponding 
to (fl m ,fl g ) = (0.242,0.703), which leads to fix = 0.055. More details on the text. 


where A = V” -SL R = V" Ei . Hi C = V” 
where A - ^ i=1 B _ ^ i=1 C _ ^ 


and Ei = -j 


The result of such analysis can be seen on Figure [tj As can be seen, the results from 
H{z) data alone yield nice constraints on the plane fl m - fl g . The flatness limit, which 
corresponds to fl m + = 1, can be seen as an straight line on this plane (dashed line on 

Fig. GO). Points on and above this line were not considered, as they correspond to non-open 
models. One may see that the best fit relies right below this line, indicating that H[z) data 
alone favour a slightly open Universe. 

Furthermore, we have considered the prior fl m > fib, with the baryon density parameter, 
fib , estimated by Planck and WMAP: fib = 0.049 0|j_a value which is in agreement with 


Big Bang Nucleosynthesis (BBN), as shown on Ref. 25]. As a result of this prior, the 3a 
c.l. contour alone is cut for low matter density parameter, as we may see on Fig. [0 

The minimum y 2 was xf nin = 16.727, yielding a y 2 per degree of freedom y 2 = 0.523. The 
best fit parameters were fl m = 0.242^g;og5lo;55 5 ^o; t | ) 9 0 , fl g = 0.703ig;34 9 lo!62 5 ia96) f° r 68.3%, 
95.4% and 99.7% c.l., respectively, in the joint analysis. 

As expected, this result is in agreement with ACDM constraints, as this model mimics the 
concordance model. Moreover, the best fit values of fl m and fl g leads to fix = 0.055, which 
corresponds to a negative value of K , as expected for this model. Sharov and Vorontsova 


6 | have found, for ACDM: fl m = 0.276^QQgg, fl\ = 0.769 ± 0.029, for la c.l., where they 


have combined H(z) with SN la and BAO data. Given the uncertainties on massive gravity 
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FIG. 2: Plot of H(z)x z for the best fit values of Massive Gravity theory (red line, from this paper, 
H(z) data only) with Tier (red dotted line). The ACDM model is also represented (black line, best 
fit from Ref. |(|, fRM+BAO+SNs). The points with error bars are the Sharov and Vorontsova 
observational data [gl. 


parameters above, the results are in good agreement, even considering the open Universe 
restriction for massive gravity, while ACDM has no restriction on curvature. 


V. BOUNDS ON THE GRAVITON MASS 


Having obtained the best fit values for the parameters, we show in Fig. [2] the plot of 
Massive Gravity theory (red line) with ±ler limit (red dotted line). The ACDM according 
to best fit data of Sharov and Vorontsova [6:] are also represented (black line). 

This model also gives an expression to the graviton mass depending on the ck 3 and ct 4 
parameters through j3±: 


^2 _ ^9 tt2 

m, - . 


(25) 


If we fix some of the parameters, we can see how the mass depends on the others. In Figure 
[3] we show a typical mass dependence with ck 4 when we set the a 3 as a constant and we 
choose to work with /3_. Such behaviour is also observed for others positive and negative 
values of a 3 . It is easy to see that the mass increases and diverges for some specific value of 
a 4 , corresponding to the limit /3_ -A 0. In some cases the mass can also abruptly decrease 
to zero, as shown in the cases cx 3 = 6, a 3 = 4 and a 3 = 2. This shows that the graviton 























mass is strongly dependent on the a’s parameters, although in the limit of very negative a 4 
values the graviton mass goes to m g ~ Hq 1 . 



FIG. 3: Some typical mass dependence with 04 for some specific values of 03 . 


VI. CONCLUSION 

In this work we have analysed an open FLRW solution of massive gravity which admits 
accelerated expansion of the universe, in full concordance to observations. The constraints 
with recent observational H(z) data were found and the best fit values obtained for the 
cosmological parameters are (fi m , Q g , Qk) = (0.242, 0.703, 0.055), in accordance with the 
ACDM model, and also point to a nearly open curvature (K < 0), for which the model is 
valid. The graviton mass dependence with the constant parameters a 3 and aq of the model 
were also analysed, and we have verified a strong dependence with such parameters. We 
have also obtained that the condition rn g ~ H ^ 1 seems dominant for a long range of the 
parameters 03 and 04 , although cosmological observations cannot be used to determine such 
parameters. 
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